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^ ' Abstract 

We study a two-dimensional fermionic QFT used to model ID strongly corre- 
lated electrons in the presence of a time-dependent impurity that drives the system 
out of equilibrium. In contrast to previous investigations, we consider a dynamic 
'^ ■ barrier switched on at a finite time. We compute the total energy density (TED) 

of the system and establish two well defined regimes in terms of the relationship 
Q ' between the frequency of the time-dependent perturbation Q and the electron en- 

ergy to. Finally, we derive a relaxation time t/j such that for times shorter than tji 
CS| . the finite-time switching process is relevant. 
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Low dimensional field theories are objects of sustained interest. In particular they are 
useful to describe the behavior of strongly anisotropic physical systems such as organic 
conductors [l], charge transfer salts [2j, quantum wires jHj, and carbon nanotubes jl|. In 
this condensed matter context, one of the most widely studied ID theory is the so called 
"g-ology" model |I31, which displays the Luttinger liquid behavior characterized by spin- 
charge separation and interaction-dependent power- law correlation functions ^. Many 
of the advances in the understanding of this model have been accomplished through the 
use of the renormalization group and the bosonization procedure, and in the context of 
equilibrium situations, i.e. in the absence of time-dependent interactions. However, as 
it is well known, in nature, thermal equilibrium is more an exception than a rule. In 
fact, in recent years there has been an increasing activity focused on the study of out of 
equilibrium low-dimensional electronic systems [7j . Thus, it becomes natural to consider 
the formulation of field-theoretical methods that can be applied to this kind of systems. 
To this aim, we have recently shown how to combine the Closed Time Path (CTP) method 
of Schwinger and Keldysh pi with a path-integral approach to bosonization ^ in order 
to analyze the g-ology model in the presence of a time-dependent dynamic barrier pTOj . 
On the other hand, as far as we know, all previous works on ID electronic systems with 
time-dependent impurities (including ref. fU]) use a large-time approximation in order 
to treat the correlation functions. In this case the effect of a non adiabatic switching 
of the interaction can be disregarded and the total energy density (TED) becomes an 
harmonic function of time which temporal average can be easily defined and computed. 
An obvious and drastic consequence is that no transient process can be predicted. The 
main purpose of this work is to reconsider this problem by taking into account the role 
of a sudden, finite-time switching process. This will enable us to show how the TED 
evolves from short to long-time regimes. 

We are interested in computing the TED for a Luttinger system defined by 

S = Sq + Sint + Sj, (1) 

where 5*0 is the unperturbed action (in the condensed matter context it is thought of as 
a linearized free dispersion relation): 



^0= / rf\^i^^, (2) 

Jc 

where J^ indicates that the integration is defined along the time contour usually defined 
in the CTP formalism [H|. Sint describes the forward scattering of spinless fermions 
(electrons): 

S^nt = -\j d'xd'y (^7/.*)(x) 1^m(x - y) (^7/.*)(y) (3) 

where x = {x, t) and the fermionic currents ^7^^? are coupled through distance-dependent 
functions, V(^)(x — y). In terms of these potentials one can make direct contact with 



the forward-scattering sector of the "g-ology" model currently used to describe different 
scattering processes characterized by coupling functions gi, g2, g^ and g^ 0. Neglecting 
processes associated to large momentum transfers, only the forward scattering couplings 
g2 and g^ play a role. The relation between these strengths and our potentials are given 
(in Fourier space, p = (p, u) ) by 



^2(P) =|(V^(1)(P) + V^(0)(P)) 
^4(P) =|(1^(0)(P)-1^(1)(P)). 



(4) 



From now on we shall disregard the momentum dependence of the potentials, i.e. only 
short-range interaction will be taken into account. The action Sj describes the interaction 
between the electrons and a localized time-dependent perturbation: 

Si= f d^x^'yoA{x,t)'^, (5) 

Jc 

where 70 is a Dirac matrix and A{x,t) is a function that contains the details of the 
perturbation, i.e. the way in which the interaction is switched on in time and the form 
in which it is localized in space. Please note that, for simplicity, only forward scattering 
between the electron and the impurity is considered. As we shall see later, this is not a 
bad approximation if weak couplings and short times are taken into account. Besides, 
in some experimental arrangements one can reasonably argue that backscattering effects 
will not be important (see, for instance ^Jj). To be specific we choose a particular form 
for the perturbation: 

A(x, t) = X sin(fit) e(t) [e(x + a/2) - e(x - a/2)] (6) 

where G(t) is the Heaviside function. We then have a separable harmonic perturbation 
suddenly switched on, consisting of a barrier of width a and height A that oscillates 
in time with frequency fl. As we shall see, this choice allows us to make contact with 
previous works on the effect of time-dependent perturbations on one-dimensional systems. 
Indeed, a similar potential (with an additional static term) was considered by Biittiker 
and Landauer in their study of the traversal time for tunneling ^2] . More recently, and 
in the context of Luttinger liquids, temporal harmonic perturbations were analyzed for 
both forward ^3] and backscattering fT^ impurities. Let us stress that neither of these 
authors includes the effect of switching, which is our main motivation. 

Let us now focus our attention on the TED, which in the Wigner representation can 
be expressed in terms of the correlation function as 

/■oo 

n{uj,x,t) = —i I (irexp(ic<jr) (^-^(x, X, t + r/2, t — r/2), (7) 



where the right and left ferinionic propagators are time-ordered along a time contour C: 
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The subscripts + and — refer to fields defined in the upper and lower branches of C, 
respectively, corresponding to forward (+) and backward (— ) time evolution. From now 
on we shall restrict our analysis to the TED corresponding to right-movers. Similar 
results can be obtained for left-moving particles. 

Using the technique described in ^U] we can factorize the Green function as 

Gc{x, y, t, t) = Gl{x -y,t~t) exp [(3{x, t) - P{y, t')] , (9) 

where the function (^^(x — y,t — t ) is the equilibrium propagator for a Luttinger liquid 
with its characteristic interaction dependent exponent 7 = (1/4) (K + K^^ — 2), {K = 

l+g4/7r-g2/7r \ ^^^ 
l+g4,/Tr+g2/Tr ' 

/3(x, t) = ^e(t){e^^*F(x, a, fi) + e-'^*F(a;, a, -fi)} (10) 

with 

F(x, a, VL) = Q{-{x + a/2)) exp[-iQ{x + a/2)] - 0(-(x - a/2)) exp[-iQ{x - a/2)] + 

+Q{x + a/2)-e{x-a/2). (11) 

Since we are mainly concerned with the effect of the switching process, we start by 
studying the free case in the presence of the time-dependent impurity (7 = 0). For this 
case the correlation function is given by 

gL°1(x, t + r/2; X, t - r/2) = ^^--. (12) 

Expanding the exponential in eq. (jH)), expressing the /3-dependent factor in terms of the 
binomial expansion and replacing in ((Zj), after some algebra one can write the TED 



as 



■-iA\" 



n{u, X, t)^=o = 0(^)+5Z {-^YY. in- ■)!■! (^(^' ''' ^))" ' (^^'^' ''' ~^))' exp{^fit(n-2j)} : 



n=l j=0 

i 



Q{uj+{n-2j)n/2)+—{Ci[{2tMao){uj+{n-2j)VL/2)]-i^i[{2t+iao){uj+{n-2j)VL/2)] + 

ZTT 



+ (-l)"+^(Ci[(2t - zao)(-cu + (n - 2j)fi/2)] - i Si[(2t - zao)(-cu + (n - 2j)n/2)])) 



+ 



cx> . , n— 1 j n—j 



n=2 j=l A:=0 /=0 ^"^ ^ ^ -^ ^ 

X ^— ^ exp{ifit(n ~2k- 21)} [Ci r(2t + iao){uj + {n - 2j -21 + 2k)n/2)] + 
27r L '- 

-iSi[{2t + iao){uj + {n-2j-2l + 2k)n/2)]-Ci[{-2t + iao)iLO+{n-2j-2l + 2k)n/2)] + 

+ i Si[(-2t + iao){uj + {n - 2j -21 + 2A;)fi/2)] 1 . (13) 

Let us mention that the first term in the above formula is the equilibrium TED. Please 
observe that for t -^ — cx), i.e. for long times in the past, only this first term remains 
and the TED recovers its equilibrium value: n{uj,x,t)^=o = no(a;) = Q{uj), as expected. 
We stress that the inclusion of the factor G(t) in the perturbation gives rise to the 
appearance of non harmonic contributions to the TED. These contributions manifest 
through the functions cosine-integral (Ci) and sine-integral (Si). For large positive times 
only harmonic terms survive. In this regime the TED is a harmonic superposition of 
equilibrium TED's centered in integer multiples of Q/2, with coefficients that depend on 
the geometry and strength of the time- dependent interaction. Only in this regime a time 
average (over the period of the interaction) does not depend on the temporal interval 
and leads to an averaged TED which is a superposition of equilibrium TED's centered 
in integer multiples of Q [13j lOj. Some of these features can be globally seen in figures 
1 and 2, where we show the behavior of n{uj,x,t)^=Q as function of uj and t, for fixed x 
inside the barrier {—a/2 < x < a/2). For comparison purposes, in figure 1 we disregarded 
the effect of the non adiabatic switching, whereas in figure 2 the B(t) function has been 
included. We have set A = 2, fi = 1 and a/2 — x = 1. 

When one focus the attention on the behavior of the TED as the energy uj is varied, 
another interesting feature is revealed. We can establish two well defined regimes in terms 
of the relationship between the external frequency Q and the electron energy u. Indeed, 
for Q/uj ^ 1 and positive times, the electrons "do not see" the harmonic perturbation 
and the corresponding TED is equal to uqIuj). So, this is the large energy region where 
the effect of the impurity is negligible. At this point we should emphasize that the 
impurity considered in this paper does not include a static term. Had we considered 
such a contribution we would have found that the large energy regime is dominated 




Figure 1: : n{uj, t)^=o as function of u and t, without finite-time non adiabatic switching. 




Figure 2: : n{uj,t)^=Q as function of uj and t, with finite-time non adiabatic switching. 




Figure 3: : An{t) in the low energy regime, for j = 0, Q = 1, u = 0.01. 



by the static impurity. The transient effects arising in this case have been extensively 
investigated in connection to the X-ray problem and the orthogonality catastrophe, for 
both Fermi ^^ and Luttinger liquids jfTO]. In our framework, however, in order to make 
a precise description of the X-ray problem, one should modify the temporal behavior of 
the impurity, i.e. the function associated to the switching process, such that a switching 
off time is also incorporated. 

Going back to the analysis of our TED, one sees that for Q/uj ^ 1 it is greatly 
affected by the impurity. In the absence of switching (figure 1) the electrons gain or lose 
energy quanta of value ±nQ (this is clearly seen when one considers the time average). 
The same occurs when switching is taken into account (figure 2), for sufficiently large 
times. Of course, there is also a crossover region given by Q/uj = 1. At this point we can 
make contact with previous studies on the traversal time for tunneling in the context of 
the quantum mechanics of a particle in the presence of a dynamic barrier jT^j. From the 
uncertainty principle we know that the traversal time r satisfies r ~ l/i^;, which leads to 
tQ = 1 , a result which is consistent with the findings of ref . J2] • 

In order to perform a quantitative description of the transient process which is now 
accesible due to the introduction of a finite-time switching mechanism, we found useful 
to consider the difference between TED's with and without the inclusion of the temporal 
Heaviside function. We call this difference An{t). In Figure 3 we plot An{t) in the low 
energy regime {Q = l,ci; = 0.01). From this figure one can see that it is possible to 
define and compute a relaxation time tn such that for times larger than tn the TED 
coincides with the one obtained when the impurity is always switched on. Thus, t < tji 



defines a temporal region in which the switching process is relevant. We have analyzed 
the envelope of An{t) for many different values of u, always in the low energy region 
Q/uj ^ 1. In other words, we fit the envelope with an exponential decay 6exp {—t/tn). 
This procedure allowed us to find t^ = l/(2uj). The variation of A and Q only affects the 
values of b. 

Let us now say a few words on the role of backscattering from the impurity. If we 
add to the Lagrangian the corresponding term, of the form ^ Abackixji)"^, of course, 
the problem of determining the Green's functions is not exactly solvable any more and 
one is forced to make a perturbative calculation. For a general interaction of the form 
Aback{x,t) = Xback f{x,t), whcrc the new coupling constant Xtack is assumed to be of the 
same order of A, it is easy to show that up to first order in these perturbative parameters, 
only the forward scattering term affects the behavior of the TED. This means that in 
the short-time regime [t < t^) and for weak coupling, only forward scattering from the 
impurity contributes to the TED. In the large-time regime, the first order contribution 
in A gives a zero average and then one is led to a higher order calculation in which 
backscattering effects would eventually contribute in a non trivial way. A quantitative 
discussion of these effects are beyond the scope of the present work, but will be the 
subject of future investigations. 

We have also computed the electrical density current J in the presence of the time 
dependent interaction studied throughout this work. Using a standard definition of J 
[T7] we obtained the following remarkably simple expression: 

J{x, t) = CX e(t) (ei-x - ^) sin[fi(t -x-^)]- e{-x + ^) sin[fi(t -x + ^)] + 

+ e{x + ^) sin[fi(t + x + ^)]- e{x - ^) sm[n{t + x- ^)]) , (14) 

where C is a renormalization constant. Let us stress that this is an exact result (in 
the absence of backscattering no perturbative expansion is required). We see that the 
current originated by the impurity (note that no external voltage has been applied) 
is a simple superposition of harmonic contributions and therefore its temporal average 
vanishes. If one introduces a bias V one obtains a (time-independent) term linear in 
V plus a photocurrent contribution similar to the above expression but with additional 
factors that depend on the right and left chemical potentials. The temporal average of 
this photocurrent is still zero. We then conclude that the dc conductance of the system 
is not affected by the perturbation. Again, this result is not changed by backscattering 
from the impurity, at least up to first order in the couplings. 

Turning back to the TED, we have extended the previous analysis to the case of a 
Luttinger liquid. Performing the same kind of manipulations as in the non interacting 
case, we found an expression similar to the one corresponding to the 7 = case (eq. (HSl))- 
Since it is another lengthy expression we shall not write it down here. The only relevant 
difference is the appearance of decaying factors of the form exp{—A/v\uj ± nQ/2\) (A is 
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Figure 4: : An{t) in the low energy regime, for 7 = 1/2, Q = 1, u = 0.01. 



an ultraviolet cutoff and v the renormalized velocity), characteristic of the smoothing 
of the energy density as consequence of the electron-electron interaction. One can also 
verify that n^jLo < "^7=0- Another interesting point we want to mention concerns the 
comparison between the areas under the curves of n and no (TED's with and without the 
time-dependent impurity) as functions of t for both the interacting and non-interacting 
electrons. In the first case the presence of the impurity leads to a smaller value for the 
total area, whereas a non vanishing forward scattering (7 7^ 0) produces the opposite 
effect. This occurs with and without a finite-time non adiabatic switching. 

In Figure 4 we show An{t) for 7 = 1/2, A = 0.1 and v = 1. We kept the same 
values used in Figure 3 for the other parameters. One observes that An{t) takes much 
lower values due to the fact that, in the presence of electron-electron forward-scattering, 
the TED's are smaller as consequence of the decaying factors mentioned above. On the 
other hand, comparing figures 3 and 4, it becomes apparent that the relaxation time is 
drastically diminished by the interaction. Then we conclude that for a Luttinger liquid 
with sufficiently large 7 the effect of the finite-time switching can be safely neglected. 



In summary, we have considered the effect of a harmonic time- dependent perturbation 
on the TED of a ID fermionic system. We put special emphasis on the case in which 
the perturbation is switched on suddenly at a finite time, and compared our results with 
the ones obtained when the dynamic interaction is present at all times. The low-energy 
transient process by which the TED corresponding to finite-time switching evolves to 
the long-time (harmonic) TED is most clearly characterized by the difference An{t). By 
carefully analyzing this function for different energies u we were able to determine that 
the relaxation time t^ associated to the TED evolution, for free fermions (7 = 0), is 
equal to l/(2a;). We also established that for Luttinger liquids (7 7^ 0) t^ is significantly 
shortened. This suggests that, in order to observe some consequence of the transient 



process, it is convenient to consider Luttinger systems with values of 7 as small as possible. 
Concerning the transport properties, we showed that the electrical current originated by 
the time-dependent perturbation is a simple superposition of harmonic contributions. 
This means that if one introduces an external voltage V (through appropriate chemical 
potentials) the dc conductance will remain unchanged by the dynamic perturbation. 
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